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OUTLINE 


I. Project Description 




ii. Turbulence Modeling 




HI. Computational Engine / Results 


FUTURE WORK 


I. PROJECT DESCRIPTION 


1- Flows of Interest 

2. Motivation 

3. Method 
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Schlieren photograph of a shock-wave turbulent boundary-layer interaction 
M=0.90 Re=1 ,750,000 [Liepmann] 


I ^MOTIVATION 





Physics 


- Boundary Layer Separation & Wall Heat Transfer 

— Spreading rate 


• Modeling 

Account for Compressibility Effects 
on Turbulence 


• Numerics 

Compare 1-point Closures on Identical Solver 


1.3. METHOD 


• 1-Point Closures: from EVM to Second-Order Closures 

• Dynamical Compressibility Effects 

• 3D / Finite Volume Approach 
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II. TURBULENCE MODELING 



1. Closure Levels 

2. Compressibility Effects 

3. Shock Wave Interactions 
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II. 1. Closure Levels 


1. EVM Mixing-Length 
(Baldwin-Lomax) 


2. EVM Multi-Equation 
(k-ff-S) 


3. Second-Order Closure 
(Shih and Lumley) 


II. 2. Compressibility Effects 
1. New Physics &c Averaging 


2. Models 
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II. 2.1. New Physics 
(Turbulent Kinetic Energy Sink) 


• e<t = (mb + jM) < d 2 > 


• nd = < pd > 


II. 2.3. Turbulence Modeling 
(Zeman, Sarkar et ai., Yoshizawa) 

• dilatation dissipation: 

£ d — (VB + ■jm) < d 2 > 

— Sarkar et al. (asymptotic analysis) 

— Zeman 

(Shocklet model) 
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• pressure-dilatation correlation: 


n d =< pd > 

- Zeman (acoustic model): 

- Sarkar et al. 

(DNS asymptotic analysis) 



Response of turbulence kinetic energy to the passage through shock 
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II. 3. Shock Wave Interactions 


1. Experimental Observations 

2. Physics 

3. Modeling 


II.3.1. Experimental Results 

• Oscillation increases with Shock Strength 
(Dolling) 

• Oscillation increases with Separation Region 

• Normal Stresses Preferentially Amplified 
(D£lery et al.) 
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II. 3. 2. Physics 


Oscillation Caused by (?): 


• ’’Breathing” of Separation Region 


• Vortex Bursting 
in Incoming Boundary Layer 
(Dolling) 


II. 3. 3. Shock Oscillation Modeling 


• Parametrized Source Terms 
in Normal RS Evolution Equation 
(gradient activated) 


• Separation region Extend 
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1. Numerical Method 

2. Turbulence Models 

3. Validation Procedure / Results 


111.1. Numerical Method 





III. 2. Turbulence Models: 

Incompressible / Compressible: an additive approach 

• Baldwin-Lomax 

• k-Epsilon / k-Epsilon-S: B.C's 

• Second-Order Closures 
Boundary Conditions: Wall-Functions 


III.3. Validation Procedure / Results 


• Calibration against simple well-documented 
flows (flat plate, jet ) 


• Results and Comparison of models 




FUTURE WORK 


Numerics 

- 2D => 3D 

- More Complex Wall Functions 

- Realizability Conditions (SOC) 

Modeling 

- Refinement of Existing Models ( e d , <pd > ) 


- Shock Oscillation Model 






Homogeneous Turbulence (R) / k-eps / Mach=0.045 Re=24357 

Downstream Evolution 
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<c <1 TKE 

v- 7 epsilon 



y=2e-5 x A (il.29) 

1 


10‘ 7 r 



v= (2e-5) x^-2.33) 
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Homogeneous Turbulence (R) / k-eps-S / Mach=0.045 Re=24357 

Downstream Evolution 


<6 <1TKE 

V — V epsilon 



Homogeneous Turbulence (R) / RSC / Mach=0.5 


Downstream Evolution 
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